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Abstract 

Strongly disordered and strongly interacting quantum critical points are difficult to 
access with conventional field theoretic methods. They are, however, both experimen¬ 
tally important and theoretically interesting. In particular, they are expected to realize 
universal incoherent transport. Such disordered quantum critical theories have recently 
been constructed holographically by deforming a CFT by marginally relevant disorder. 
In this paper we find additional disordered fixed points via relevant disordered deforma¬ 
tions of a holographic CFT. Using recently developed methods in holographic transport, 
we characterize the thermal conductivity in both sets of theories in 1+1 dimensions. 
The thermal conductivity is found to tend to a constant at low temperatures in one 
class of fixed points, and to scale as T 0 3 in the other. Furthermore, in all cases the 
thermal conductivity exhibits discrete scale invariance, with logarithmic in temperature 
oscillations superimposed on the low temperature scaling behavior. At no point do we 
use the replica trick. 
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1 Introduction 


1.1 Universal incoherent thermal transport 

The thermal conductivity is the simplest and most universally defined quantity in which 
the subtleties of dc transport arise and can be studied. In particular, in a translationally 
invariant system at a nonzero temperature, the thermal conductivity is generically infinite. 
This is because at a nonzero temperature there will be a nonzero energy density and hence an 
overlap between the thermal current and the momentum operator. Momentum is conserved 
by assumption and does not relax. Hence the thermal current also cannot fully relax. 

The fate of momentum conservation is especially important in strongly interacting sys¬ 
tems. In a weakly interacting theory, there are many long lived quantities, the quasiparticle 
number excitations rik- Momentum is a particular linear combination of the rq c and as such 
is often not especially privileged. However, in the absence of quasiparticles, but in the pres¬ 
ence of translation invariance, momentum conservation dominates the late time behavior of 
heat current excitations. 

Strongly interacting transport in the presence of weak momentum relaxation has been 
extensively studied in recent years. In particular, the holographic correspondence has pro¬ 
vided tractable explicit models of strongly interacting dynamics. One important upshot 
of this work has been the development of the old memory matrix formalism mm as the 
framework of choice to discuss non-quasiparticle transport with weak momentum relaxation 
MEUamElEUinittl]. In particular, if the low energy physics is described by an IR fixed 
point with an emergent long wavelength translation invariance (as is the case whenever the 
low energy dynamics admits an effective gapless QFT description), then the memory matrix 
tells us that the dc thermal conductivity is controlled by the leading irrelevant operator that 
breaks translation invariance [a si- Transport in these cases is now solved in principle, 
up to the characterization of the leading irrelevant operator for a given system. 

A more challenging scenario is a strongly interacting theory in which translation in¬ 
variance is not approximately restored even at the lowest energy scales. In this case the 
memory matrix is not a useful tool. If the system does not become thermally insulating 
then the low temperature thermal conductivity is an intrinsic property of the low energy 
physics. In the context of charge transport, such systems were called ‘universal incoherent 
metals’ in [T2j. It was suggested in [12] that the cuprates and other ‘bad metals’ should 
be understood as such universal incoherent metals. The objective of the present work is 
to provide a controlled theoretical model of universal incoherent thermal transport. The 
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central role played by disorder means that the model is unlikely to be directly relevant to 
the cuprates, although it may have features in common with superfluid-insulator or quan¬ 
tum hall plateaux transitions. The objective is to gain a theoretical and conceptual handle 
on incoherent transport. Other recent work towards a theoretical framework for incoherent 
transport, especially charge transport, can be found in mm- 

1.2 Disordered quantum critical points 

Known instances of gapless low energy theories without translation invariance are disordered 
quantum critical points (or phases). These arise when disorder is a relevant deformation 
of a quantum critical theory mm which does not gap the theory, but rather the system 
flows to a new critical point with a finite amount of disorder (yet another possibility are 
so-called ‘infinite disorder’ fixed points HZmEUHUB], these will not be considered here). 

In principle, controlled interacting disordered fixed points could be obtained by applying 
epsilon or large N expansions directly to a replicated held theoretic description of a known 
CFT, such as the Wilson-Fisher fixed point. However, it is found that the beta function 
equations do not admit perturbative, stable zeros. Instead, one finds a runaway flow to 
strong disorder mm- The only known quantum critical disordered fixed points found 
by this method involve either continuing to a small number of time dimensions [2Tj or 
allowing disorder with long range correlation in time [22] • It is unclear if these additional 
approximations can fully capture the quantum disordered dynamics as they are essentially 
expansions close to the classical statistical physics disordered fixed point [23]. In any case, 
since these fixed points are found in perturbation theory, transport will be described by a 
quasiparticle-based Boltzmann equation, and will not access the strongly interacting and 
strongly disordered regimes we are interested in. 

Recent work has used the holographic correspondence [23J to obtain controlled, strongly 
interacting and strongly disordered fixed points |25, 26]. Aspects of these fixed points will 
be recalled below. They were obtained by deforming a CFT by marginally relevant disorder. 
Past work showed that the (line of) fixed points are characterized by a dynamical critical 
exponent z > 1 that becomes larger as the tunable disorder strength V is increased. In 
parallel with these developments, major recent work has essentially ‘solved’ the problem 
of computing d.c. conductivities in holographic theories [27J [28j, 29 ] EDJ. In particular, 
the new results give a powerful way to compute the thermal conductivity of complicated 
disordered spacetimes without needing to explicitly solve perturbation equations about the 
background. In this work we obtain the thermal conductivity of the holographic disordered 
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fixed points found in [25,[26]. We will also construct new disordered horizons arising from 
relevant rather than marginally relevant disorder, and obtain the thermal conductivity of 
these new fixed points. 


1.3 Results 

The main results in this paper will be for disordered fixed points in 1+1 dimensions. There 
is no obstruction to considering higher dimensions beyond the fact that more computing 
power is needed to solve nonlinear partial differential equations in more dimensions. The 
results for the thermal conductivity are as follows: 

• At weak disorder, in both the marginally relevant and relevant cases, the thermal 
conductivity k(T) is excellently described by the perturbative formulae ([5]) and 
over the entire temperature range that we can access. See figures [l] and [5] below. 

• At stronger disorder, the thermal conductivity sees the onset of a low temperature 
scaling regime. This regime is characterized by a discrete scale invariance, so that: 



( 1 ) 


For exponents a and j3, a scale To and a function F(x) « co + c\x, for the regimes 
considered. See figures [2] and [6] below. 

• For the cases of marginally relevant disorder, our results are consistent with the scaling 


exponent in © taking the value a = 0, independently of the dynamical critical 


exponent z in the low energy fixed point theory. The magnitude Co of the conductivity 
decreases with increasing strength of disorder. For the relevant disorder, we find 
a « 0.3. This latter exponent appears to be independent of the strength of the 
disorder in the UV and hence a property of a disordered fixed point theory. 

To our knowledge, these are the first theoretically controlled results on incoherent transport 
in an explicit strongly disordered and strongly interacting system. In addition 

• For the newly constructed disordered fixed points from relevant disorder we obtain 
the entropy density as a function of temperature. This thermodynamic quantity is 
also found to exhibit discrete scale invariance at low temperatures. See figure |4j 

As anticipated, the thermal conductivity at low temperatures is given by a power of 
temperature that is intrinsic to the disordered quantum critical low energy theory. These 
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systems therefore indeed exhibit universal incoherent thermal transport. Interestingly, how¬ 
ever, the power of temperature that appears - a in ([Tj) - is not the exponent (14) below 
that would be anticipated from the simplest dimensional analysis. This is presumably pos¬ 
sible because of the scale To in the low energy theory (or perhaps, due to a broad disorder 
distribution at the disordered fixed point), and suggests that scaling analyses for incoherent 
systems, such as m, should proceed with caution. 

The exponent a = 0 in the marginally relevant case indicates that the conductivity 
takes a finite universal value in these cases (up to oscillations in the logarithm of the 
temperature). The value depends on where we are in the line of disordered fixed points, 
but is an intrinsic property of the low energy theory. A universal thermal conductivity is 
reminiscent of [32], where the strength of the disorder appears in two quantities (the finite 
quasiparticle lifetime and the finite density of states) that cancel in the thermal conductivity. 
That mechanism, however, is rather perturbative. It seems more likely that the universality 
here is tied to having d = 1 spatial dimensions and the fact that the fixed point is obtained 
as a marginally relevant deformation of a z = 1 UV theory (from which, perturbatively, 
a temperature independent thermal conductivity is obtained ([8])). By this logic, in d = 2 
spatial dimensions we would expect to obtain k ~ T in the marginal cases. 

Discrete scale invariance has appeared in several previous descriptions of disordered 
fixed points pi m ssi m m, in the running of couplings towards the fixed point theory. 
The appearance here in a radically distinct and nonperturbative theoretical framework - 
and, crucially, without use of the replica trick - suggests that this is indeed a robust feature 
of disordered fixed points. Complex scaling exponents in holography are often indicative of 
dynamical instabilities [36], or that the quantum phase transition is pre-empted by a first 
order transition m- However, the aforementioned perturbative results show that stable 
spiraling or limit cycle renormalization group flows are possible in disordered systems. 


2 Disordered horizons and disordered quantum criticality 

In this section we summarize the holographic description of disordered fixed points | 25l [26] . 
In doing so, we also recall several general concepts of disordered quantum criticality. 

The starting point will be a conformal field theory (CFT) in d spatial dimensions. This 
theory is deformed by an operator coupled to a quenched random potential 

£ = £cft + V(x)0(t, x) . (2) 

The operator O has mass scaling dimension A while the disorder is taken to be short range 
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and Gaussian, so that the disorder averages 

(V(x)) R = 0, (V(x)V(y)) R = V 2 d^ ( x - y ). ( 3 ) 

The Harris criterion mm states that the coupling ([2]) is then relevant if A < (d + 2)/2. 

A minimal holographic setup that can describe a CFT deformed by the random coupling 
<§ is gravity coupled to a scalar field, with action 

5 = 16 ^ / ^ [ R + " 2 (V< ^ )2 “ 4 C/($) ) ' (4) 

Following the usual holographic dictionary [24] . a potential of the form [/(T) = —/jl $ 2 /(2L 2 ) 
means that the bulk scalar field <f> is dual to an operator O with dimension A satisfying 
A(A — d — 1) = —/t. In the works [25., [26] we set /i = d(d + 2)/4, so that A = (d + 2)/2 
and the disorder is marginal, saturating the Harris criterion. In particular, this means that 
the disorder strength V in ^ is a dimensionless free parameter. In section [6] of this paper 
we will also consider the case of relevant disorder with A = 3/4ind=l dimensions (for 
context, marginal disorder has A = 3/2 in d = 1, and the unitarity boundary is A = 0). 

The papers |25l [26] asked: what is the low energy dynamics of the deformed theory (§? 
In holography this question is studied [24] by solving the Einstein equations of motion follow¬ 
ing from the action Q subject to the boundary condition that the scalar field tend to the dis¬ 
order potential V at the asymptotic boundary of the spacetime: lim r _^o [Y A_d_1 cI>(r, x)] = 
V (x). The Einstein equations with a disordered source are complicated. They can be solved 
numerically or perturbatively in weak disorder. With the solution at hand, the low energy 
physics - such as entropy density at low temperatures - is characterized by the properties 
of the spacetime deep in the interior, where the holographic coordinate r —> oo. The re¬ 
sults of [25., [26] will be reviewed below. Other works constructing disordered holographic 
spacetimes include [Ml [39], 001 HU 02] 03]. 

The disordered potential V(x) is modeled by a sum over N cosines with increasing 
wavevector and random phases. We will not give details here, see [25i f26]0 There is both 
a short and long distance cutoff on the wavevectors, that we will denote as k\jy and &ir 
respectively. Note that, /cir = k\jy/N. This translates into a temperature range of roughly 
&ir < T < fcuv over which the system is truly disordered. Thus, accessing disordered 
physics at low temperatures requires a large number N of oscillatory modes. In fact, as 

X A small difference in notation compared to the earlier papers: in this paper we use d to denote the 
number of boundary spatial dimensions. Also, we have corrected a factor of 27r in the normalization of V 
in j25j [261 , so that now V 2 = V 2 OIlect = 2-7rV/j d . To facilitate comparison with previous works, especially for 
numerical results, we will sometimes quote the result in terms of V = Hid, so that V 2 = 27rT 2 . 
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explained in [26], one is able to go to slightly lower temperatures than the rough estimate 
suggests thanks to favorable numerical factors. 


3 Perturbative formula for the thermal conductivity 


In a regime where the disorder may be treated perturbatively, analytic held theoretic results 
are possible using the memory matrix formalism. For a recent clean discussion of the 
hydrodynamics of a strongly interacting CFT perturbed by weak disorder see [ 41 ]. The 
thermal conductivity is 


s 

K pert. = p i 


( 5 ) 


where s is the entropy density, and the momentum relaxation rate mum 

r-Ylr f ddk fc 2 Im G$ 0 (u,k) 
sT J (2n) d d uj 


( 6 ) 


Here d is the number of spatial dimensions and Gqq(uo, k ) is the retarded Green’s function 
of the operator that couples to disorder in ([2]). For the case of strongly interacting theories 
without an underlying relativistic scale invariance see mm- 

From the above two formulae one can immediately obtain the temperature scaling of 
the disorder at temperatures &ir < T < fcuv (this is roughly the range of temperatures that 
is insensitive to the cutoffs on the disorder distribution). Simple power counting gives [4] [5] 

j'2d+l-2A 

K pert. ~ ^- ■ (7) 


Here we used the fact that s ~ T d for a CFT. In particular, for the marginal case 


ftpert. ~ 


1 

~w~ 


( 8 ) 


These expressions use the scaling symmetries of the high energy CFT fixed point. They need 
not hold at the lowest temperatures, even for small V, as there the system has been driven 
to the disordered IR fixed point by the relevant or marginally relevant deformation ([2]). 
Understanding the thermal conductivity in the far IR is the primary objective of this paper, 
and we will return to that question shortly. At weak disorder, the perturbative thermal 
conductivity ([T]) can be expected to appear at intermediate temperatures. Reproducing the 
perturbative answer will build confidence in our more general results. 

Specializing to the holographic theory allows us to compute the retarded Green’s func¬ 
tion Gqq(uj, k ) in ([g]) explicitly. In this computation and in the remainder of the paper we 
will focus on one boundary space dimension (d = 1). Both the holographic formula for dc 
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conductivities and the numerics we shall employ below are substantially simpler in this case. 
There is no fundamental obstacle, beyond computing power, to study of the higher dimen¬ 
sional case. The details of the (somewhat standard) holographic computation of Gq 0 (uj, k) 
are given in Appendix [A} The perturbative momentum relaxation rate is found to be 


T = 


V 2 L 1 

2ttGn sT 


(27 tT) 


2A _ 3 sin 2 (7rA) 


7 


r(2 - A) 2 


f kvv dk , 2 


f A 

ik \ 

/ — k 2 

T I 


+ 

Am 27r 


V 2 

4t tT 


(9) 


We have explicitly included the long and short wavelength cutoffs on the disorder in order for 
better comparison with the numerical results below. The perturbative thermal conductivity 
([5]) is then found by recalling that for the three dimensional Schwarzschild-AdS black hole 
that we are perturbing about, the entropy density s = (ttL)/(2Gn) x T. 


4 Holographic formula for the thermal conductivity 

In a sequence of beautiful recent works, Donos and Gauntlett have obtained a general non- 
perturbative formula for the thermal and electric dc conductivities of holographic theories 
mmmm- Their expressions are a major generalization of the older important result 
by Iqbal and Liu [45j. Specifically, for a class of holographic theories, they have derived an 
expression for the dc conductivities in terms of data evaluated purely on the event horizon. 
Crucially for our present purposes, the formula holds with arbitrary spatially dependent 
sources in the dual quantum field theory. These results allow us to go beyond perturbation 
theory in the disorder. 

In this section we adapt the results of (25] to the case of the Einstein-scalar theory 0 in 
three bulk spacetime dimensions. Scalar fields have also recently been considered in [Ml ED]. 
The details of the computation can be found in Appendix [B| The full spacetime metric is 
written as 

L 2 r 


ds z = 


y 


X 2 \ , B(y,x) 2 

+ FW dy 


( 10 ) 


y 2 + y—f(y)A(y,x)dt 2 + S(y,x) (dx + F(y,x)dyY 

Here f(y) = 1 — y 2 . This form of the metric will be useful for the numerics in the following 
section. In these coordinates, the horizon is at y = 1 and the asymptotic boundary at y = 0. 
At the boundary we impose A = B = S = 1 and F = 0, so that the only source will be 
in the scalar field ~ r d+1 ~ A V(x). At the horizon A = B = A^°\x), S = S^°\x), the 
scalar $ = and F again vanishes. With these boundary conditions, the constant 

y+ determines the temperature as y+ = 2nT. In Appendix [B] the thermal conductivity is 
found to be 


k = 


7 t 2 T 2 

2Gn 


1 

k T 


dx 


(tw^r 

L\/sm 


-1 


( 11 ) 
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Here L x is the length of the x direction in the dual field theory. 

To evaluate the integral in ( |11[ ), we must solve the Einstein-scalar equations of motion 
numerically in order to find the disordered black hole background and hence <h®(x) and 
S'®(x). The virtue of (11) is that we do not need to additionally solve numerical pertur¬ 
bation equations about the background. In Appendix [B| we check that ( |TT| ) reproduces the 
perturbative memory matrix result ([5]) at weak disorder. Schematically: To leading order in 
small V, S'® = 1 and describes a linearized perturbation on top of Schwarzschild-AdS' 3 . 
The integral f dx (9 X $®) 2 appearing in (11) can then be written in terms of the low fre¬ 
quency limit of the retarded Green’s function of O appearing in ([ 5 ]). In fact, it is of note 
that the full holographic expression © looks like a local version of the perturbative result 
([ 5 ]), in which S'® plays the role of a spatially varying entropy density. This ‘hydrostatic’ 
aspect of holographic dc conductivities has recently been emphasized in [T3j . 


5 Numerical results 

Consider first the case of marginal disorder: A = 3/2 in d = 1. For this case, we have 
constructed the numerical solutions using the same methods as we described in [26]. The 
spacetime takes the form ©• We solve the Einstein-scalar equations following from the 
action Q, subject to the boundary conditions described below equation ( |To] ) above. From 
the solution, the thermal conductivity is obtained from horizon data using ©• All numerics 
in this section are performed with a sum over N = 50 oscillatory modes in the disorder 
potential V{x). 


5.1 Thermodynamics and scaling expectation for the conductivity 


It is found that in the CFT deformed by marginal disorder, the disorder grows logarith¬ 
mically towards low energies (i.e. towards the interior of the spacetime), and is therefore 
marginally relevant [25l [26] . Such logarithms were first found holographically in [39] . 

A basic quantity characterizing any disordered fixed point is the dynamical critical 
exponent z. In [25] it was found that the disorder averaged spacetime metric in the far 
interior, at zero temperature, takes a ‘Lifshitz ’ m scaling form 


lim (ds 2 )n ~ 

r—¥ 00 



dr 2 + dx 2 
r 2 


( 12 ) 


The dynamic critical exponent z was determined analytically in a perturbation theory in V 
and numerically in general. Furthermore, heating up by a small temperature T, the entropy 
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density was shown to scale as 


s ~ T d/z , (13) 

with the same z as appeared in the disorder averaged metric @0 These last two results 
above demonstrate the emergence of a disordered quantum critical point at low energies, 
with an emergent scale invariance characterized by the exponent 2 . The fact that z > 1 
indicates that disorder is playing a finite role in the fixed point physics. 

Before turning to numerical results for the thermal conductivity, the thermodynamic 
facts above suggest the following simple scaling analysis. The entropy scaling s T d / Z 
suggests that hyperscaling is obeyed and that there are no scales remaining at the disordered 
IR fixed point. In particular, this scaling suggests that the energy density e has its canonical 
dimension [e] = d + 2 . The conservation equation e + V ■ = 0 then determines the 

dimension of the heat current to be \j®\ = d — 1 + 2 z. The thermal conductivity is defined 
via jQ = —acVT, and hence [k] = d + z — 2. We conclude that 


Scaling ~ T^ d+Z ~ 2 ^ Z T 1 - 1 / 2 


(14) 


This result gives k ~ T° when z = 1, consistent with the perturbatively marginal disorder 
expression (|8j) when d = 1. However, for any arbitrarily small disorder strength we have just 


recalled that the IR fixed point has z > 1. The scaling result (14) would therefore suggest 


that the thermal conductivity should go to zero like (14) at the lowest temperatures, as 
the marginally relevant disorder drives the system away from the UV fixed point. This 
expectation will be seen to fail. 


5.2 Thermal conductivity 

Figure [l] shows the thermal conductivity as a function of temperature for weak disorder, 
V = V/\/2tt = 0.01. Also shown in the same plot is the perturbative result ^ and 
The numerical and perturbative results are seen to be in excellent agreement over the entire 
temperature range. We have made a log log plot to emphasize the two different power law 
regimes. The high temperature power law asymptotes to k ~ T 3 at temperatures T > k\jy, 
i.e. above the short distance cutoff on the disorder distribution, as can be seen from the 

2 It should be noted that the difference between the scaling form s ~ T x l z (in d = 1) and the perturbative 
result s ~ T (1 + [4 — l] logT + ■ ■ ■) can be quite small when z ~ 1. In particular, in perturbation theory, 
the logarithmic derivative ^^f; = (; + ( 3 — l) 2 logT + • • • is close to the exponentiated value of 4, even at 
low T. However, at V = 1 one has z ~ 2 [2S], and the hence the constant logarithmic derivative found in 

1261 requires the exponentiated scaling behavior of the entropy density. 
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perturbative result ([9]). In the truly disordered regime, the thermal conductivity is seen 
to be a constant. This regime extends over 0.01 < T/k\jy < 0.06 in the plot. There is 
no sign over this temperature range of the IR cutoff on the disorder distribution. While 
&ir = k\jy/N = 0.02 k\jy here, numerical factors shift the effect of the cutoffs to lower 
temperatures [26] . 



^UV 

Figure 1: Log-log plot of thermal conductivity k(T) with V = V/\[2t: = 0.01 and 
marginal disorder (A = |, d = 1). The dots are numerical data points with N = 50 oscillator 
modes generating the disordered potential. The solid line is the analytic perturbative result 
([ 5 ]) and ([ 9 ]). At low temperatures k ~ T°. At high temperatures compared to the disorder 
cutoff k\jY, k ~ T 3 . 

Bolstered by the agreement of the previous figure, we can now turn to stronger disorder. 
Figure [ 2 ] shows the thermal conductivity as a function of temperature for V = 0.5 and 
V = 1. The plots focus on the low temperature regime, away from the short and long 
distance cutoffs on the disorder. These are the temperatures over which the scaling of the 
entropy density s ~ T 1 / 2 was found in [2SJ • Once again, the thermal conductivity saturates 
to a constant value at low temperatures. The value of the low temperature conductivity is 
no longer close to the perturbative prediction. A new phenomenon appears, furthermore, 
which is that the thermal conductivity oscillates at low temperatures. 

In the following section we will consider relevant disorder. For the relevant disorder 
we will be able to go to lower temperatures and see more oscillations than are visible 
in figure [2j The perseverance of the oscillations with stronger disorder suggests that the 
oscillations continue down towards zero temperature and are not merely an intermediate 
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Figure 2: Low temperature thermal conductivity k(T) for marginal disorder with 
V = 0.5 (left) and V = 1 (right). The disorder has been simulated with N = 50 oscillatory 
modes. The perturbative formulae ([5]) and Q predict the zero temperature values of the 
conductivity to be 0.72 and 0.18, respectively, and therefore underestimate the true result. 

scale phenomenon. The simplest explanation is that they correspond to an imaginary scaling 
exponent for the thermal conductivity, which leads to log-periodic oscillations (discrete 
scale invariance) of the form advertised in equation ([!]) above. We mentioned previously in 


section 1.3 above that discrete scale invariance has appeared in previous perturbative field 


theoretic studies of disordered fixed points, but had been suspected to be an artifact of the 
replica trick. We are not using any replica trick here. We will do a detailed fitting to the 
parameters in 0 in the following section on relevant disorder, where we are able to get to 
lower temperatures. 

The fact that the conductivity tends to a constant (up to oscillations) at low tempera¬ 
tures is in disagreement with the simplest scaling expectation ©■ For instance, equation 
(14) would predict k ~ T 0 ' 5 for the V = 1 case, for which z ~ 1 


This is not 


seen in figure [2] Instead the results in figure [2] suggest that the low temperature scaling 
of k ~ T° is correctly predicted by the perturbative formula ([8]), up to the appearance of 
oscillations and a renormalization of the overall value of the conductivity. At present we do 
not understand the physics at work here. 


6 Relevant disorder 

This section gives evidence for a new class of disordered fixed points, obtained by deform¬ 
ing a CFT by relevant rather than marginal disorder. Because these are new gravitational 
solutions, we will describe how they have been constructed and characterize their thermo¬ 
dynamics before going on to present the thermal conductivity. 


12 






6.1 Numerical methods 


For sufficiently negative mass squared, there are two choices of quantization in asymptoti¬ 
cally Anti-de Sitter spacetime. Motivated by the expectation that disorder effects should be 
stronger for the choices of A that most violate the Harris criterion, we decided to work in 
alternative quantization. For reasons we explain shortly, we chose /x = 15/16 in the action 
Q which, using alternative quantization, leads to A = 3/4. This in turn means that, in 
Fefferman-Graham coordinates, the scalar field approaches the boundary (z —> 0) as 

$>{z,x) = z 3//4 [{0{x)) + yfzV(x) 4-], (15) 


where V(x) is the disorder potential. 

We will be using spectral collocation methods which, in their simplest formulation, are 
best tailored for solving PDEs where all variables are analytic in their domain of integration. 


This suggests the line element (10) should be changed in order to accommodate the slow 
decay exhibited by the scalar field in (Il5|). We choose: 


ds 2 =£ 


y+ \-g(y)A(y,x)dt + S(y,x) (dx + 2yF(y,x)dy)‘ 


4y 2 B(y,x) 2 

+-tvt- dy 

9{y) 


(16) 


where g(y) = 1 — y 4 . If we set A = B = S = 1 and F = 0, we recover the BTZ black hole. 
For the scalar field we set <3?(y, x) = y 3 / 2 y?(y, x). In order to solve the resulting system 
of PDEs, we follow mutatis mutandis [26] . In particular, we use the De-Turck trick [48] , 
with the BTZ black hole being the reference metric. The only change worth mentioning 
is the boundary condition for the scalar field (or alternatively, for ip). We first note 
that, asymptotically, the relation between y and the usual Fefferman-Graham coordinates 
is simply y = yd/I/y^fl + 0(z 3 / 2 )]. This implies that if we want to recover the expansion 


(15) we need to set a Neumann-type boundary condition for ip of the form 

dy(p(y,x)\- 0 = y~ b/i V{x). 


(17) 


The variables are analytic in these new coordinates. The fact that analyticity can be 
achieved through a relatively simple change of coordinates is the reason we focus on A = 3/4. 

In the numerics for this relevant case, the disorder potential will be modeled with N = 50 
or 100 cosines. In performing the actual numerics one can set ko = 1, which effectively sets 
the units. Since the Einstein-scalar equations are nonlinear, we want to make sure we can 
capture the first five harmonics emanating from each progenitor k n mode. In order to make 
sure we resolved all such scales, most of our simulations run with 100 points in the radial 
direction y and no less than 1000 points in the Fourier direction x. An illustrative solution 
is shown in figure [3] 
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Figure 3: Numerical solution with relevant disorder, V = 1 and N = 100. The 

plots show the disorder potential V(x), the induced expectation value (0(x)) for the scalar 
field, the scalar field on the horizon and the scalar field everywhere in the spacetime. 

6.2 Thermodynamics 

This subsection will describe the dependence of the entropy density s on temperature in two 
cases with relevant disorder (A = |) and with a fairly strong disorder magnitude in the UV: 
V = V/\/2tt = 0.5 and V = 1. In the case of marginal disorder, logarithmic divergences 
in perturbation theory could be plausibly resummed. Relevant disorder leads to power law 
divergences in perturbation theory that should not be naively resummed. Therefore we 
must turn to numerics from the outset to access the low temperature thermodynamics. In 
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appendix [C] we characterize, in the spirit of [25], the perturbative power law divergences 
that appear in the averaged zero temperature metric towards the interior of the spacetime. 

A plot of entropy versus temperature in these theories shows an essentially linear in 
temperature entropy over all temperatures. This is the result for a non-disordered one 
dimensional CFT, and suggests that z remains equal or close to one. However, a more 
detailed probe at low temperatures reveals more dramatic structure. Figure [4] shows the 
logarithmic derivative of the entropy with respect to temperature at low temperatures. This 
is the same diagnostic that was used in [26] to conclude that s ~ T 1 /^ at low temperatures 
in the marginal case. Figure [4] shows that instead, in this relevant case, the numerical data 
is well fit by the discrete scale invariant form 


T ds 
~sdT 


jn 


^b 0 + b\ sin 



(18) 


The values of the parameters appearing in the fit are given in the figure caption. The 
specific heat, c = T ds/dT, is always positive in these plots, as required for thermodynamic 
stability. We do not quote To, as it is ambiguous up to To —> e 27 rn/,< 5 To. The precise 
form of this function is not the most important issue. Rather, what the fits in figure [4] 



Figure 4: Logarithmic derivative of the low temperature entropy density s(T) for 
relevant disorder with V = 0.5 (left) and V = 1 (right). The disorder has been simulated 
with N = 100 oscillatory modes. Dots are numerical data whereas the solid line is a fit 
to the discrete scale invariant form (18). The values in the fit are {7 ~ 0.02 ,6 
1 , b\ 


—0.02} for V = 0.5 and then {7 « 0.03, 5 ~ 4.2, 6 q ~ 1, 61 


4.5,6o 

—0.04} for V = 1. 


show is that there is log-oscillatory structure in the thermodynamics of these models. The 


frequency of the oscillations, 6 in (18), will be seen to be close to the frequency (3 of the 


log-oscillations found in the thermal conductivity in the following section. The small value 


of 7 in (18) together with 1 ~ 6 q 3> &i is behind the fact that the entropy is almost linear 
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in temperature over this range of temperatures. The nonzero 7 exponent suggests that a 
different temperature scaling may set in at much lower temperatures. The fact that the 
parameters 7, <5, bo are similar in the two cases supports the hypothesis that the solutions 
are flowing towards the same IR fixed point. 

It seems possible that such oscillations are also present in the entropy density of the 
marginal case at sufficiently low temperatures. 

6.3 Thermal conductivity 

We start with a weakly disordered case to corroborate our expectations. Figure [5] shows the 
thermal conductivity as a function of temperature for weak disorder, V = V /\/2vr = 0.01. 
Also shown in the same plot is the perturbative result ^ and (J9|^J now with A = |. Once 
again, the numerical and perturbative results are seen to be in excellent agreement over 
the entire temperature range. The log-log plot cleanly shows the three expected behaviors. 



^UV 

Figure 5: Log-log plot of thermal conductivity k(T) with V = V/\[2t: = 0.01 and 
relevant disorder (A = |, d, = 1). The dots are numerical data points with N = 50 oscillator 
modes generating the disordered potential. The solid line is the analytic perturbative result 
([ 5 ]) and ([ 9 ]). At high temperatures compared to the disorder cutoff k\jy, K T 9 / 2 . In 
the intermediate, disordered regime k ~ T 3 / 2 . At the lowest temperatures, below the long 
wavelength cutoff fciR, n increases, ultimately exponentially. 

At temperatures above k\jy, the thermal conductivity tends to n ~ T 9 / 2 , as can be seen 

‘To capture the very low temperature regime correctly, i.e. T ~ km, one must replace the integral in ([ 9 ]) 
with a discretized sum over N terms. 
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from the analytic expressions. At temperatures below fciR, the potential is effectively a 
lattice with wavevector &ir rather than a disordered potential. A lattice cannot efficiently 
relax momentum in a relativistic theory and so the conductivity is expected to increase 
exponentially towards low temperatures J6jj. The intermediate truly disordered temperatures 
show the scaling k ~ T 3 / 2 anticipated from perturbation theory ([t]) . 

Figure [6] shows the thermal conductivity for stronger disorder, V = 0.5 and V = 1. The 
plot focusses on the disordered low temperature regime, but above the IR cutoff (these plots 
have used N = 100 oscillatory modes to simulate the disordered potential, allowing us to go 
to lower temperatures). Also shown in the plot is a fit to the discrete scale invariant form 


^d.s.i. 


rjiCX 


^c 0 + ci sin 



(19) 


The values of the parameters appearing in the fit are given in the figure caption. The 
numerical data is seen to contain over one period of oscillation, close to two periods for the 
V = 1 case. There are some error bars on the fits, due to the finite temperature range, 



^UV 


Figure 6: Low temperature thermal conductivity /c(T) for relevant disorder with 
V = 0.5 (upper curve) and V = 1 (lower curve). The disorder has been simulated with 
N = 100 oscillatory modes. The solid lines are fits to the discrete scale invariant form of 
equation (19). The fit for the V = 0.5 case has {a ~ 0.3, /3 ~ 5.4, cq « 3,ci « 0.06}. The 
fit for the V = 1 case has {a ~ 0.28, /3 « 4.1, cq ~ 2, ci ~ —0.09}. 


but both disorder strengths are consistent with roughly a ~ 0.3. Note that this exponent 
is significantly different to the perturbative exponent of 3/2 that we saw in figure [5j This 
scaling is therefore indicative of a new strongly disordered fixed point. Also in both fits 
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ci -C Co, consistent, perhaps, with the truncation to a single oscillatory frequency. The 
agreement in the values of the exponent (3 is less good - this is likely because of the limited 
temperature range of the oscillations in the V = 0.5 fit. 

We see that although the magnitudes of the thermal conductivity are comparable to 
those for marginal disorder in figure [2]- over this range of low temperatures - the tempera¬ 
ture dependence is quite different. The conductivity in the relevant case is tending to zero 
at low temperatures whereas the marginal case appears to be stabilized. Note that we do 
not expect any general lower bound on the thermal conductivity analogous to that found 
in mm for the electrical conductivity. This is because the thermal conductivity does not 
have a ‘pair production’ term. A natural quantity to consider in incoherent transport is the 
diffusivity [12] 

D = - K — . 

Tds/dT 

For the cases we have considered, the diffusivity indeed grows both towards low and high 
temperatures, with a minimum at intermediate temperatures. 

7 Towards possible applications 

Disorder is likely to be an important actor in many real world quantum critical systems. 
Some of the most interesting putative quantum critical physics in condensed matter systems 
- for instance in the cuprates or heavy fermions - involves a nonzero charge density (i.e. 
‘metallic’ quantum criticality) and is beyond the scope of the framework we have consid¬ 
ered. Three important instances of ‘zero density’ quantum criticality - in which in addition 
disorder appears to play a key role - are (i) Superfluid-insulator transitions mm, (ii) 
Metal-insulator transitions [51, 152] . (iii) Quantum hall plateaux transitions [53], 154] . The 
thermal conductivity is not straightforward to measure in these systems, but would offer 
direct insight into the role of disorder in the quantum critical physics. 

The systems described in the previous paragraph have two spatial dimensions. While 
in principle the methods we have used adapt to two space dimensions, the additional com¬ 
puter power required, especially to get down to low temperatures with strong disorder, is 
substantial. The real challenge here is to obtain as explicit an analytical handle as possible 
on disordered horizons, beyond perturbation theory. Disordered zero temperature horizons 
may exist as decoupled fixed point solutions of the Einstein-scalar system in the same way 
that extremal horizons describe IR fixed points. Understanding such solutions on their own 
terms would amount to substantial progress in the theory of disordered quantum criticality. 
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A Green’s function of the scalar operator 


This appendix computes the low frequency Green’s function for a scalar operator in the 
background of a neutral black hole in AdS 3 . We can write the background as |24] 

ds2 = ~ r 2 (-/( r )^ 2 + J ^) + dx 2 ) ’ ( 21 ) 

with / = 1 — r 2 /r^_. This spacetime is a solution to the theory Q with no scalar turned 
on. The temperature is T = l/(2nr + ) and the entropy density s = (ttL)/(2Gn) x T. 

The equation of motion for the scalar field $ about this background is 


V 2 d> = 


A( A — 2) 
L 2 


d>. 


( 22 ) 


We have written the equation in terms of the dimension A of the dual operator O. The 
slick argument in | 10 ] gives the imaginary part of the retarded Green’s function at low 


frequencies in terms of a non-normalizable solution of (22) with no time dependence. That 
is, write = $o(u k)e lkx and then [10] 


lim 

CJ—>-0 


Im Gq 0 (uj , k ) 


T oo 

oj 


L 1 


$ 0 (r+,k) 2 . 


(23) 


47 tGat r+ 

Here <I>o(r, k ) must be taken to be regular at the horizon and to have the non-normalizable 
near-boundary behavior of < I , o(r, k) ~ r 2 ~ A . The prefactor of 47 tGtv is related to the 
normalization of the scalar action Q . For equation (22), the relevant solution is 


Mr,k) = r 2 ~ A 2 F 1 I 1- 


A 2—2A r(2 — A) 

r 7+ r(A) 


ik 


ik 


A 

4t tT ~ ~2 ,1 + 4t tT 


- ^,2-A,(2vrTr ) 2 


( 24 ) 


r (A 4 . jai 

1 I 2 ^ 47 tT) 


r (!-f +&) 


„ , ik A ik A . . 9 

2Fl 1 "UrT + ¥’4^T + Y ,A,(27rTr) 
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Therefore from (23) 


lim Im G% 0 {u,k) 

a;->0 OJ 


L 


AttGn 


(2 vrT) 


2 A— 3 sin 2 (tt A) 


7T 


r(2 - A) 2 



(25) 


B Horizon formula for the dc conductivity 


This appendix adapts the computation in [28] to obtain the thermal conductivity of Einstein- 
scalar theory in three bulk spacetime dimensions. 

As described in the main text, the background metric takes the form 


2 r 


ds 2 = — 2 ~ 

y 


V+ ( ~f(y)A(y, x)dt + S{y,x) (dx + F(y,x)dy) ) + 


B(y,x)dy 

f(y ) 


2 "I 


(26) 


with / = 1 — y 1 . The scalar field takes the form <J> = #(y, x). In the main text we described 
the boundary conditions satisfied by the various functions in the metric. To obtain the 
thermal conductivity one can consider perturbations of this background of the form 


ds 2 —> ds 2 + h a b(y, x)dx a dx b + 2t(g tt dtdx . 


(27) 


The f-linear source in 5gt x is equivalent to a thermal gradient d x T = —T ( 1273 - 


In the perturbed metric (27), construct the tensor 

G^ = V [/ ^ ] , k = 


8 ( 

dt 


(28) 


Note that k is not a Killing vector in the perturbed spacetime. The manipulations in [28] , 
with the additional requirement that the scalar field be independent of t, show that the 
quantity 

Q = V-gG xy , (29) 


is independent of both x and y. Using the metric (26) one obtains 

A 3 / 2 


Q = 


i/l_ 

L 2 y B 1 / 2 S 1 / 2 


d v 


y 2 h tx 

fA 


- 8x 


y 2 h ty 

fA 


(30) 


From this expression it can be verified that d x Q = d y Q = 0 using the equations of motion. 

The fact that Q is constant throughout space means that it can be evaluated both on the 
horizon and at the boundary. Equating the two resulting expressions will give the thermal 
conductivity. Near the boundary one finds 


UmQ = 8irG N (T tx ) FT + 0(t). (31) 

y-> o 

Here (T tx )p t is the heat current (also, the momentum density) in the dual field theory, 
computed as usual from the holographic stress tensor [55! ]. Note that to from the bulk 
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stress tensor T ab to the field theoretic momentum density (T tx )p^, one has to multiply the 
bulk T ab by y/a£ a nb, where n is a unit timelike vector, £ = d x and the spatial measure 


a = y/ g xx ■ The time dependent part in (31) is unimportant for our purposes f271, 28]. 
Now we consider the near horizon behavior of the fluctuations. We only need the 


behavior of the two metric components that appear in the expression (30) for Q. The 
general behavior of these components near the horizon is 

h tx {y,x) = Ly + S (0 \x) 1 h[°J(x) + h[l\x)f(y) log {jXXy) + °( y ~ ^ ’ ( 32 ) 


ht y {y , x ) = hfj(x) + 0{y - 1 ). 


(33) 


The logarithmic term is included in ht x so that the total 5gt x = tQgu + Sht x can be written 
in terms of the ingoing Eddington-Finkelstein coordinate v, which is: 


v = t - 1 log 
2 y+ 


i +y 
i -y 


(34) 


Regularity (infalling boundary conditions) requires the t and y dependence of 5gt x to com¬ 
bine into a dependence only on v. Therefore 


hS(x) = L 2 A {0 \x). 


(35) 


Plugging the near horizon behavior into Q, we find Q = —y + h[^(x) + 0(y — 1). Since 
this must be a constant, we see that d x iif^ = 0. We can go further and set all of the higher 
order terms (in y — 1) to zero, i.e. impose d y Q = 0. The leading term in this equation near 
the horizon reads, using the background equations of motion expanded near the horizon: 


0 = y+C + 


2 h® (d x $wy 


- d x 


h 


(o) 

ty 


h 


(o) 

tx 


- r d x log A® 


(36) 


Ly + Ss$) \_L 2 AW 2 Ly+VS^ 

If we average this equation over a period L x in the x direction, the total derivative term 
drops out. This allows us to solve for h!fj. Insertion into the expression for Q gives 

,3 


d >( o) (u 

Q = - y+Kx = 


2X 


Since y+ = 2vrT, the thermal conductivity is: 


L x J L\/l 5 ^) 


K = 


(T 


'xt\ 


/FT 


1 


Q 


IT 


■T 2 1 


-d x T 8v tG n CT 2 G N X 


(37) 


(38) 


The second equality used the expression (31) for Q at the boundary while the third equality 


used the expression (37) for Q at the horizon. Equation (38) is the result for the thermal 


conductivity quoted in the main text. 
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Let us check the perturbative limit. Here S ^ = 1 and is the horizon value of 
the solution of a linearized scalar about the Schwarzschild-AciSs solution with the random 
source asymptotic boundary condition [25' [26]: 

N 


Ak 

<b(r, x) = 2VJ — ^2 ^o(r, k n ) cos (k n x + y n ), 


2vr 


(39) 


n =1 


where A k = k\jy/N, k n = nAk, the phase is uniformly randomly distributed between 0 
and 27 t, and the solutions at each k n are normalized so that ‘bo ~ r 2 ~ A near the boundary, 
with unit coefficient. Therefore 

r2 N 


X = T - 


kj f ’ (&4>) 2 = L£(Ai)^4. 0 (r + ,fc 

Jx n=l 


(40) 


Now using (23) and taking the continuum A k —> 0 limit 

x = 4' G N y 2 r°° dJ ^ k 2 lim Im Gq 0 (uj, k) _ nT ^ 
L 2 T J 0 2tt uj->o lo L 


Here T is defined in ([b]). Hence (38) becomes 


s 

K= f 


(41) 


(42) 


Where we used s = ( ttL)/(2Gn ) x T. We have reproduced the memory matrix result ([5]). 
To get the match we have included a factor of 27r in the normalization of (39) that was 
missing in 


C Perturbative backreaction of relevant disorder 


In this appendix we consider the lowest order backreaction of relevant disorder, with A < 
(d + 2)/2, on the metric. We work in d = 1 and at zero temperature. The calculation 
proceeds very similarly to that in (25) . so details will be sparse. The final result is equation 


(48), which shows that two components of the averaged metric have power law growths as 


one moves towards the interior of the spacetime. The power law growth is of course as one 
expects for a relevant operator. A minor novelty compared to the marginal case (25] is that 
two rather than one metric component diverge (in Fefferman-Graham gauge) towards the 
interior of the spacetime. 

The scalar solution that corresponds to a disordered source on the boundary was given 


in (39). To lowest order at zero temperature, the scalar propagates in an AdS 3 background. 
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Therefore the profile with the required near-boundary behavior is a linear combination of 
Bessel functions: 


23-A / Afc 

<f>(r, x) = V V y, rk£~ 1 K A - 1 (k n r ) cos ( k n x + 7 n ) 


(43) 


where k n , 7 „ and AA; are defined below (39). Relevant disorder for d = 1 corresponds to 
A < 3/2. 


The scalar field (43) sources metric perturbations at order U , through the Einstein 


equations following from the action 0 - The metric to second order can be written: 

1 


ds = ^ 


(44) 


— ^1 + V 2 A^ 2 \r, x)j df 2 + dr 2 + ^1 + V 2 B^ 2 \r, x)J dx 2 

A/ 2 ) and are fixed, modulo boundary conditions, in terms of the scalar sources via the 
Einstein equation. The solution reads as: 

A^ 2 \r,x) = ai(r) + a 2 (x) — 4 J drdx d r <&d x A , (45) 

B (2 \r, x) = /?i(x) + ^/3 2 (x)r 2 + 2 J dr r j ^ [/J?4> 2 - r 2 (d x 4>) 2 - r 2 (<9 r <f>) 2 ] , (46) 

1 ( f dr 

Q?i(r) = r/i + - r] 2 r 2 + 2 drr / (^ 2 4> 2 + r 2 ^^) 2 - ?’ 2 (<9 r <f>) 2 ) . (47) 

Regularity at the Poincare horizon r —> 00 requires r] 2 = a 2 = /3 2 = 0, while j3\ is fixed by 

setting A ^ = B^ at the conformal boundary r —> 0. 


Upon plugging in the scalar solution (43) and averaging over disorder configurations, 
i.e. the random phases 7 n , one can readily extract the large r behavior of the averaged 
metric components from these integral expressions. We find: 


lim (A^Ur, x)\ =rji—'-/V 2 r 3 2A , lim (B < ' 2 \r,x)\ = /3i(x) — 6V 2 r 3 

■—>00 \ / R r—»• 00 \ / R 

r(A-i)r(2A-§) 


—2A 


where 7 and <5 are the positive constants 

(3 - 10A + 4A 2 )T (A - §) T (2A - f) 


7 = 


6 = 


4 A T(A — 1) 2 T(A + 1) ’ 4 A_ 1 T(A — 1) 2 T(A) 

3 


(48) 


(49) 


Note that as A —> the marginal limit, 7 has a pole, and one recovers the logarithmic 


divergence found and resummed in [25] to obtain a non-trivial dynamic critical exponent. 
Away from this limit, the power law divergences cannot be resummed reliably. 
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